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A continuum-space bead-spring model is used to study the phase behavior of binary blends of 
homopolymers and the structure of the interface between the two immiscible phases. The structure 
of the interface is investigated as a function of immiscibility, chain length, and system size. Capillary 
waves are observed and their measurement allows us to determine the surface tension 7. We propose 
a more universal method of measuring the interfacial width in terms of second moments of the 
different contributions to the first derivative of the interfacial profile. Predictions of this method 
are directly verified. The effect of chain length on the surface tension is also studied. 

PACS numbers: 61.41+e, 64.60. Cn, 61.25.Hq 
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The mixing of two different homopolymer species, say 
A and B, in the molten state often results in a system 
of two immiscible phases. The structure of the inter- 
face between the phases has been the object of several 
theoretical studies p]|| which use as a starting point the 
Flory-Huggins free energy of the blend and its well-known 
associated interaction parameter xab |§ ■ These theories 
predict that the intrinsic order parameter profile of an 
interface located at z a is well-approximated by [H 



ip(z) = -0 o tanh[2(z - z )/w ] 



(1) 



where ip(z) = (p A (z) - pb(z))/(pa(z) + pb{z)), pi is the 
number density of species /, and tp is the bulk value 
of the order parameter. According to the same theo- 
ries 0, the intrinsic interfacial width w is predicted 
to follow w = a(6xAs) -1 ^ 2 , while the surface tension 
7 = (ksT '/ 'a 2 )(xAB I '6) 1 / 2 . Here, a is the statistical seg- 
ment length, ks is the Boltzmann constant, and T is the 
temperature. 

On the experimental side, neutron reflectivity experi- 
ments Q have measured a value of w that is systemat- 
ically and substantially larger than that predicted theo- 
retically from the xab associated to these blends. More 
recent experiments [|| seem to support the proposed ex- 
planation H that the presence of capillary waves at 
the interface could be responsible for this discrepancy. 

Of similar importance is the effect of chain length N 
on the interfacial width, since theoretical treatments are 
mostly in the limit of infinite N . Some finite-chain- length 
corrections have been proposed ^ but none of these have 
been tested systematically. 

With these motivations in mind, a simple continuous- 
space (CS) model has been developed || in order to study 
binary blends of homopolymers. It is the purpose of this 
Letter to present results, as derived from this model, for 
capillary-wave and chain-length effects at the interface 
of immiscible homopolymer blends. Our approach differs 
significantly from previous numerical studies which use 
lattice models almost exclusively Jul. For systems with 



interfaces, CS models have interesting features: besides 
their inherent spatial isotropy and the absence of pin- 
ning, they offer a simple way to determine the surface 
tension from the measured pressure tensor [ fj"l"|] . More- 
over, provided the forces are short-ranged, theoretical 
work |l2| suggests that interfaces in the continuum ex- 
hibit no roughening transition. 

The CS model used here has been described in more 
detail elsewhere In this model, the polymer chains 
are represented by attaching N soft spherical beads of 
mass m using a finitely extensible spring potential. The 
softness of the beads, or mers, is set by the interaction 
potential. A binary system is built by constructing a 
large number M of such chains, of given type A or B, 
and enclosing them in a virtual box having the desired 
boundary conditions. 

We use molecular dynamics (MD) as the simulation al- 
gorithm, supplemented by occasional Monte Carlo (MC) 
type exchanges in order to improve the sampling of phase 
space. The time evolution of the coordinates of all chains 
is resolved through Newton's equation of motion (EOM), 
integrated using a velocity- Verlet algorithm [ JL3| . The 
motion of the chains is coupled to a heat bath act- 
ing through a weak stochastic force W(i) and a corre- 
sponding viscous damping term with friction coefficient 
r. Besides improving the diffusion of the system in phase 
space, this coupling has the practical advantage of sta- 
bilizing the numerical calculation, especially in the pres- 
ence of the energy fluctuations induced by the MC ex- 
changes. Including these terms, the EOM of mer i reads 



d 2 r % 



dt 



W,(t). 



(2) 



The last term Wi(t) is a white noise having an aver- 
age strength determined by temperature and the friction 
coefficient through the fluctuation-dissipation theorem. 
In the viscous drag term, F has to be carefully chosen 
to avoid overdamping so that the motion of the mers 
be dominated by their inertia. Finally, the conservative 
force term derives from a potential energy U having two 
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contributions: the interaction potential Uu acting be- 
tween all mers, responsible for excluded volume effects, 
and an attractive potential holding adjacent mers along 
the chain U . For the first contribution in a mixture 
of two homopolymer species, the interaction potential 
Uu{rij) between beads i and j of types I, J = {A, B} 
separated by a distance is taken as the repulsive core 
of a central-force Lennard- Jones (LJ) 6:12 potential, 



Uu(nj) = 4e/j 



12 



(3) 



for TYj < r c = 2 1 / 6 cr and zero otherwise. Instead of using 
the Lorentz-Berthelot mixing rule, we choose a = <ta — 
as, and ejj = (1 + SjjS)e, where S > is a small param- 
eter controlling the miscibility and Sjj is the Kronecker 
delta. Our choice in energy parameters can be seen as a 
special (symmetric) case of S — [e^s — \ {^aa + ess)]/e, 
which represents the kind of interactions used in simple 
lattice systems such as in the Flory-Huggins theory. Here 
e and a are, respectively, parameters fixing the energy 
and length scales. Accordingly, our results are reported 
in terms of these natural units, using r = a{m/e) 1 / 2 for 
the time scale. All the results reported here are obtained 
using r = 0.5T" 1 at a constant overall number density 
p = 0.85cr~ 3 (isochoric) and a temperature T = I.Oe/fcs. 
Chain size N varies from monomers to 30-mers. 

Simply by varying S, while leaving the temperature 
constant, one can induce a phase transition from a ho- 
mogeneous blend to a system with two coexisting phases. 
This effect is reminiscent to what is done in recent exper- 
iments reported by Gehlsen et al. ]l4| ] in which the phase 
separation of binary mixtures is studied as a function of 
the difference in deuterium between two otherwise iden- 
tical polymer species. Thus, in principle, the model and 
an MD simulation method would be adequate to study 
the kinetics of phase separation in binary blends. How- 
ever, the time scale involved in the unmixing transition 
— occurring through the diffusive motion of the chains — 
is extremely long. To properly sample the phase space 
in a reasonable calculation time, we supplement the MD 
moves with an MC procedure relabeling the type of the 
homopolymer chains (A <-» B) with a Metropolis tran- 
sition rule |Xq] . Our exchange attempt rate is approxi- 
mately M chains per r M. 

Cubic simulation cells of volume L 3 = NM / p and with 
full periodic boundary conditions were used to study the 
bulk properties of our model 0. The equilibrium bulk 
order parameter ip a (M, TV, 5) was measured for different 
sizes as a function of the immiscibility parameter 5. A 
finite-size scaling analysis Q of our results, involving the 
multiple- histogram technique [ jig ], was used to extract 
the critical immiscibility value S C (N) associated with a 
system of infinite size (M — > oo). For N > 2, our data 
is well-described by NS C = 3.40(5)fcBT/e, clearly evi- 
dencing a I/TV behavior of the critical immiscibility, as 




FIG. 1. The interfacial profiles 'i(z) defined as the xy-slice 
averages of ip(r) for a system of 4000 10-mers with antiperi- 
odic boundary conditions. The values of S c /8 are as indi- 
cated. The lines are from a least-squares fit to the function 
= Y? i= i ^erii^z/wi). 



predicted by the Flory-Huggins theory. 

The finite-size effects on the interfacial properties of 
multiphase systems are not as well-known as for the 
bulk. Capillary- wave effects 0, confinement [jl7| or 
finite-chain-size Q corrections to properties such as the 
interfacial width or the surface tension are still under ac- 
tive investigation. We shall now address these questions. 

Once the critical miscibility parameter is known, we 
can build antiperiodic immiscible systems (at S c /5 < 1) 
and study the interfacial properties of the model, that 
is, we use simulation cells having an antiperiodic [^8 |l9 
boundary condition in one direction and a volume L 2 L±, 
where L± is the length perpendicular to the interface. 
Figure |f| shows the interfacial profile of the order param- 
eter Vt'(z) of such systems, obtained by spatially averag- 
ing ip(r) in xy slices, perpendicular to the antiperiodic 
direction z, and then averaging the resulting profiles in 
phase space with the proper offset (since the zero of the 
interface is translationally invariant). It is tempting to 
fit the order parameter profile directly to the same intrin- 
sic interfacial function ^f(z) = tp a tanh(2z/w), and thus 
associate a width w to the interface. However, our aver- 
aging of the profile incorporates the effect of the intrin- 
sic interfacial width w as well as any other perturbation 
broadening the interface, such as capillary waves. To sep- 
arate both effects, a systematic study of the interfacial 
profile for different system sizes has to be performed. 

In capillary wave theory [Q, the roughness ((Azo) 2 ) 1 / 2 
of the zero of a sharp interface z (x, y) is found to be, in 
the linearized regime of small distortions, 



((Az ) 
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In 



(4) 



where is the lower cut-off driven by the system size 
and generally taken to be n/L while qi is the upper cut- 
off, usually assumed to be driven by some correlation 
length such as 7r/(c'u> ), where c' is some number. It 
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is assumed here that L± — > oo so that the roughness is 
strictly controlled by L. 

We shall now derive a general and formal argument al- 
lowing us to add the roughening effects of capillary waves 
to the intrinsic interfacial width, which we shall measure 
universally as a moment. If one assumes that capillary 
waves can be decoupled from density fluctuations, the 
averaged interfacial profile $f(z) can be written as the 
convolution of the intrinsic interfacial profile tp(z — z a ), 
and the probability V(z )dz of finding the interface at 



(5) 



(z) = / i/>(z- z )V{z ) dz D . 



By applying d/dz on each side, one finds that #'(z) is 
the convolution of two well-bounded functions, ip'(z — z a ) 
and V(z ). We associate a functional measure v[f] to a 
well-bounded function / as pO] 



v[f] 



Too/W dz 



dk 



-f(k)\ 



k=0 



/(0) 



(6) 



where / is the Fourier transform of /. For real /, it 
can be shown pU using the convolution theorem that 
v[^f'] = v[ip'] + v[V], which we rewrite as 



A = A; 
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where A 2 = v[^'] now measures the total interfacial 
width while A 2 = v[ip'} is related to the intrinsic interfa- 
cial width w a . The last term comes from Eq. |] and the 
fact that v[P] = ((Az ) 2 ). We reiterate here that Eq. 
is of general validity since it does not rely on any specific 
functional form and only uses well-controlled approxima- 
tions. For example, if one fits ^(z) to cither a simple 
tpo tanh[2z/w] or to Y^i=i cueifiy/nz/wi) |22|, then A 2 is 



2 w 2 /48 or a * w i/( 27r Z)™=i a i)> respectively. 

In Fig. we present results from fitting ^(z) to n=2 
error functions. The dependence of A 2 on both the lin- 
ear size L of the system (while L± ^> A) and 5 clearly 
evidences the presence of capillary waves. To our knowl- 
edge, this is the first time that capillary waves have been 
observed in computer simulations of polymer interfaces 
p9| . We fit our data to the expression A 2 = as + bs ln(L), 
where b$ ~ 1/7l (subscript L indicates that the surface 
tension has been derived from capillary- wave effects) . By 
measuring the surface tension "fp independently from the 
difference in the pressures (obtained from the virial jl3) , 
hence the subscript P) perpendicular and parallel to the 
interface fujl , jp = (Pj_ — P»)Lx, one gets another esti- 
mate of the surface tension, and can compare it unam- 
biguously with 7i [H). Figure || shows the two different 
measures of 7 for a monomeric system. For all N, our 
data show that with j L — k B T/(2irbs), both measures 
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FIG. 2. Linear-size and miscibility-parameter dependence 
of the interfacial width for various systems of 10-mers. Values 
of S c /5 are as indicated. The lines are least-squares fits to 
A 2 — as + bs ln(L). Systems are of size L 2 L±, with L± S> A. 
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FIG. 3. Two different estimates of the surface tension 
for monomeric systems. Surface tension yp is obtained from 
L±(P± — P||) while 7_l = kpT '/(2irbg), where bs is obtained 
from a plot similar to Fig. [| 



are consistent, thus evidencing the origin of the broaden- 
ing as being capillary waves. 

The other parameters as obtained from the least- 
squares fits are always negative. Once we know as and 
bs, Eq. |^ can be solved numerically for A D . It turns out 
that our data set only has a solution for c > 13, thus im- 
posing a bound on cut-off qi . Figure ^ shows values of A 
extracted by assuming c — 13 and plotted such as to com- 
pare with theory which predicts that w%/N ~ 1/(Nxab), 
where N\ab is the measure of immiscibility in the Flory- 
Huggins theory (with the critical value (Nxab)c = 2), 
such that 2/{Nxab) corresponds to S c /S for our model. 

The Lj_-dependence of the interfacial width has also 
been investigated. For L± > 4 A (and 5 c /5 < 0.8), we 
find no effect of L± on A, in contrast with recent studies 
]l7t that observe that A ~ L 1 / 2 for near-critical immis- 
cibilities. In these studies however, L± is such that the 
interface is confined so that the tails of the order param- 
eter do not relax to the bulk value. In the present work, 
this effect is taken care of by only using results from sys- 
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FIG. 4. The scaled intrinsic interfacial width A^/N as 
estimated from Eq. ([?]) and c = 13.. The chain length varies 
from monomers to 30-mers. The data are plotted so that the 
reduced miscibility S c /S is equivalent for all chain lengths. 
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FIG. 5. The miscibility-parameter dependence of the sur- 
face tension for systems of different chain lengths. The scaling 
form holds very well over the range studied. 



tems for which L± is large enough so that there exists 
a wide region (going through the antiperiodic boundary) 
in which ip is fully relaxed at ±ip {8). 

In Fig. ||, the effect of the chain length on the surface 
tension is shown. From self-consistent field results (and 
as N — > oo), N 1 / 2 ^ ~ (Nxab) 1 / 2 ■ Accordingly, our 
result collapse very well when N 1 ' 2 "/ is plotted in terms 
of 5 c /5. 

We have shown that a continuous-space model with 
a rather simple potential allows us to study the inter- 
face between immiscible binary blends of homopolymers. 
The natural isotropy of an off-lattice model coupled with 
a weak mismatch between the polymer chains made it 
possible for us to observe strong capillary-wave effects at 
the interface of moderately small systems and evaluate 
the surface tension unambiguously. 
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